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ABSTRACT 

We present measurements of the auto- and cross-frequency correlation power spectra of the cosmic 
(sub)milhmctcr background at: 250, 350, and 500 \mv (1200, 860, and 600 GHz) from observations 
made with the Balloon-borne Large Aperture Submillimeter Telescope, BLAST; and at 1380 and 
2030 |J.m (218 and 148 GHz) from observations made with the Atacama Cosmology Telescope, ACT. 
The overlapping observations cover 8.6 deg^ in an area relatively free of Galactic dust near the south 
ecliptic pole (SEP). The ACT bands are sensitive to radiation from the CMB, the Sunyaev-Zcl'dovich 
(SZ) effect from galaxy clusters, and to emission by radio and dusty star- forming galaxies (DSFGs), 
while the dominant contribution to the BLAST bands is from DSFGs. We confirm and extend the 
BLAST analysis of clustering with an independent pipeline, and also detect correlations between the 
ACT and BLAST maps at over 25(7 significance, which we interpret as a detection of the DSFGs in 
the ACT maps. In addition to a Poisson component in the cross-frequency power spectra, we detect a 
clustered signal at Aa, and using a model for the DSFG evolution and number counts, we successfully 
fit all our spectra with a linear clustering model and a bias that depends only on redshift and not on 
scale. Finally, the data are compared to, and generally agree with, phenomenological models for the 
DSFG population. This study represents a first of its kind, and demonstrates the constraining power 
of the cross-frequency correlation technique to constrain models for the DSFGs. Similar analyses with 
more data will impose tight constraints on future models. 

Subject headings: cosmology: cosmic microwave background, cosmology: cosmology: observations, 
submillimeter: galaxies ~ infrared: galaxies ~ galaxies: evolution - (cosmology:) 
large-scale structure of universe 
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1. INTRODUCTION 

Roughly half of all the light in the extragalactic sky 
which originated from stars appears as a nearly uniforni 
cosmic infrared b ackground (GIB; iPuget et al.l 119961 : 
iFixsen et al.lll998l ). This background peaks in intensity 
at around 200 |xm ()Dole et al.l 120061 ) , and results from 
thermal re-radiation of optical and UV starlight by dust 
grains, meaning that half of all the light emitted by stars 
is hidden by a veil of dust. 

Following its discovery, stacking analyses have statis- 
tically resolved most of the GIB shortward of 500 |xm 
into discrete, dusty star-forming galaxies (DSFG), and 
to a lesser extent radio galaxie s , at z < 3 (e.g. . 
Dole et al.ll2006HDevhn et al.l l2009l: iMarsden et 311120091 : 



Pascale et al.ll2069l ). Longward of 500 [a.m, the contribu- 
tion from radio galaxies and higher-redshift DSFG to the 
GIB hicrcascs dramatically with increasing wavelength 
(e.g.. iBcthermin et al.|[2010l ). and as a result, the GIB at 
these wavelengths has yet to be fully resolved into dis- 
crete sources (e.g.. IZemcov et aI1l2010[ ). At wavelengths 
longward of ^ 1 mm, while both radio sources and DS- 
FGs (e.g.. IWeifi erall [20091 : IVieira et al.|[2010l ) are stiU 
present, signal from the cosmic microwave background 
(GMB) becomes visible and dominates the power on an- 
gular scales larger than ^ 7 arcmin [l ^ 3000) at A = 2 
mm. 

To fully realize the cosmological information encoded 
in the GMB power spectrum, contributions to it from 
sources must be removed. At current mm-wave detec- 
tion and resolution levels, the radio sources are primarily 
discrete (Poisson) while the DSFGs arc confusion limited 
and clustered. For example, at 148 GHz the power spec- 
trum of DSFGs roughly equals the GMB power spectrum 
at ^ » 3000. Thus knowledge of DSFGs is important for 
understanding the scalar spectral index of the primordial 
fluctuations and other parameters encoded in the high-£ 
GMB power spectrum. 

Because GMB maps contain signal from multiple con- 
tributors, determining precisely the level at which galax- 
ies contribute to the GMB power spectra is non-trivial. 
Submillimeter (submm) maps, on the other hand, for 
the most part contain signal from dusty galaxies, so 
that cross-frequency correlations of submm and mm- 
wave maps provide a unique way to isolate the contri- 
bution of DSFGs to the GMB maps. However, submm 
maps of adequate area and depth have until now not ex- 
isted. 

Here we present the first measurement of the cross- 
frequency power spectra of submm and mm-wave maps. 
We use mm-wa ve data from the Atacama Gosmolo gy 
Telescope (AGT: lFowler et al.ll2007l : lSwetz et aLll2010t ) at 
1380 and 2030 fim (218 and 148 GHz), collected during 
the 2008 observing season, and submm wave data from 
the Balloon-bo rne Large Aperture Submillimeter Te le- 
scopc (BLAST : iPascalc ct al.ll2008HDevlin et al.ll2009D at 
250, 350, and 500 ^m (1200, 860, and 600 GHz), which 
were collected during its 11 day flight, at ^ 40 km al- 
titude, in Antarctica in 2006. We use these to mea- 
sure the power from DSFGs, both Poisson and clustered. 
These results will complement those anticipated from 
the Planck mission (jTauber et al.|[2010t ) by extending to 
higher resolution in the mm-wave regime. 

This paper is organized as follows: In § [2] we briefly 



overview the sources of signal in the submm and mm- 
wave sky, their spectral signatures, and the models we 
adopt to describe them. In § [3] and |4] we describe the 
data and detail the techniques used to measure the power 
spectra. We present our results in §[5l and interpret them 
in terms of a linear clustering model in § [6l We discuss 
and conclude in § [7] and [H] 

2. THE (SUB)MILLIMETER BACKGROUND 

The dominant contribution to the cosmic submillime- 
ter and millimeter-wave background, referred to here- 
after as the GSB, depends strongly on wavelength and 
angular scale. One map may have contributions from 
galaxies, GMB, and the Sunyaev-Zcl'dovich (SZ) effect 
simultaneously. 

2.1. Power Spectra 

The beam-corrected power spectrum of the sky is a 
superposition, depending on wavelength, of the following 
terms: 
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where Ci represents the angular power spectrum in mul- 
tipole space, £, and Ni is the noise. Here "cirrus" refers 
to emission from Galactic dust (§ 12. 6p . "ff" refers to 
free-free emission, and "radio" refers to radio sources, 
whose flux increases at longer wavelengths (§ 12. 4p . It is 
assumed that diffuse synchrotron emission is negligible. 
For the purposes of this paper we consider GHz-peaked 
sources and similar objects as radio sources. Equation [T] 
assumes that the various components are uncorrelated, 
when in reality, correlations among various components 
likely exist. Typically, however, these correlations should 
be small and can be reasonably neglected. Further- 
more, we define the power from the extragalactic sky as 
CP^ = Cf^ - Cf - Cf. The Cf component is neg- 
ligible for the area of the sky we are dealing with and we 
ignore it. In what follows we report cross power spectra 
as both Ce and P{kg), with ke the angular wavcnum- 
ber. To convert from multipole £ to kg, or from /iK^ to 
Jy^ sr~^, see Appendix \K\ 

In order to isolate the spectra of one or more contrib- 
utors to the background requires removal, or adequate 
modeling, of the unwanted power. Since the contributors 
have distinct spectral signatures (i.e., their flux densities 
vary from band to band differently), multi- frequency ob- 
servations make decomposition of the signal possible. For 
discrete sources, the ratio of flux densities from band to 
band is 
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inhere 



is the "spectral index" , and is a function 



of the rest-frame spectral energy distributions (SEDs) 
of the sources that make up the galaxy population, and 
their redshift distributions. Consequently, measurements 
of the spectral indices can place powerful constraints 
on source population models (e.g.. IMarsden et al.l 120101 : 
iBethermin et al.ll2010l ). 

2.2. CMB 
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At wavelengths longer than ~ 1 mm [v — 
350 GHz) the CMB dominates the power spec- 
trum on scales greater than ^ 8'. Multiple 
peaks in the spectrum ha. ve been measured most 
recent l y by: iBrown et al.1 (I 2009f): iFriedman et al.l 
(l2009l): iReichardt et all (l2009bllaD : ISavers ct all (120091): 
iLuekeret all (I2010D: iSharp et all (I2010D: iFowler et al I 



( 2O1O0 : iDas et all ()2010[ ): and iNoha et all (j2009f ) . See- 
ondary anisotropics include: CMB lensing, which acts to 
smooth out the peaks and add excess power to the damp- 
ing tail; and the SZ effect, which distorts the primordial 
CMB signal. 

In the present analysis, the CMB power in ACT maps, 
which dominates on large scales, does not correlate with 
signal in the BLAST maps; however, it does act to in- 
crease the noise on those scales (see Appendix |B]) . 

2.3. Dusty Star-Forming Galaxies 

DSFGs, as their name implies, are galaxies undergo- 
ing vigorous star formation, much of which is optically 
obscured by dust. They have average flux densities of 
5 mJy (at 250 |J^m: iMarsden et al.l [20091 ) . star-formation 
rates (SFRs) of ^ 100 - 200 Moyr^i (jPascale et all 
120091: iMoncelsi etHI 120101), number densitie s of - 2 x 
lO"'^ Mpc"^ fe.g.. Ivan Dokkum et aIll2009D . and typ- 
ically lie at redshifts - 4, with the pea k in the 
distribution at z 2 (jAmblard et al.l I2010D . They 
are distinguished from "submillimeter galaxies" (SMGs) 
discovered by SCUBA_ (|Smail et all 119971 : iHughes et all 

fi99a 

dant 



Bales et al.lll999D. which are ten times less abun- 
Coppin et al.ll2006l ). lie at slightly higher redshifts 
(jChapman et al.l [20051) . have SFR - 1000 Mgyr"!, and 
are thought to be triggered largely by major mergers 
([Engel et al.ll2010D . They are of course related: SMGs 
comprise the extreme, high-rcdshift end of the DSFG 
population. 

The dust in DSFGs absorbs starlight and re-emits it 
in the IR/submm, with a spectral energy distribution 
(SED) phenomcnologically well approximated by a mod- 
ified blackbody, 

cx v^B{i^), (3) 

where B^iy) is the Planck function, and /3 is the emis- 
sivity index (whose value typically spans 1.5 - 2; e.g., 
iDraine fc Le"a[l98l . The SED of a typical DSFG with 
temperature 30 K and (3^2 (jChapin et al.l 120101) 
peaks at rest-frame A ~ 100 i^im (which redshifts into the 
submillimeter at z ~ 1 — 10). A property of this shape 
is that with increasing redshifts, observations in the 
(sub)mm bands continue to sample at a rest-frame wave- 
length close to the peak of the SED, so that even though 
sources become more distant, their apparent flux remains 
roughly constant. This so-called negative K-correction 
makes observations at lon ger wavelengths m ore sensitive 
to higher redshift sources (jBlain et al.il2003D . As a result, 
sources at z ^ 1 have a significant impact on the power 

spectrum at (sub)mm bands. 

The cross-power spectrum between frequency bands 1 
and 2 arising from unclustered sources is related to the 
number counts, i.e., the surface density as a function 
of flux density (S) as follows 



where dN/ dS AS is the number of sources per unit solid 
angle in a flux bin of width AS, at frequency bands 1 
and 2, and 5'cut is the flux density at which the counts 
arc truncated. When the slope of the counts is steeper 
than —3, the power diverges at low fiux densities, and the 
power spectrum is dominated by the contribution to the 
background from faint sources. In the case of DSFGs, the 
strong evolution of the source cou nts with redshift res ults 
in a steep slope at the faint end ([Devlin et al.ll2009t ). so 
that after masking local sources (with S'cut ~ 500 mJy) 

the DSFG component of the GIB at A > 250 |J.m is dom- 
inated by faint sources and remains finite. 

Since galaxies are spatially correlated (being biased 
tracers of the underlying dark matter field), the DSFG 
power spectrum has both Poisson and clustered compo- 
nents: 
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The measured strength of the clustered component 
is such that it domin a tes o ver the Poisson on 
scales > 3' ([Viero et all I2009L hereafter V09) and 

(IMarsden ct all 120091 iHall et all 120101 : iDunklev et all 
l2010HShirokoff' et al.ll2010l ). withthe exact value depend- 
ing on frequency and flux cut. How the strengths of the 
Poisson and clustering terms scale with wavelength, and 
if they evolve together or independently, remain an open 
questions. 

We compare to t he models of [ Bcthermin ct al.l (|2010l 
hereafter BIO) and IMarsden et al.l ([20101 hereafter MIO) 
for DSFGs at the BLAST and ACT wavebands. These 
are phenomenological models which are specifically tai- 
lored to constrain the evolution of the rest-frame far- 
infrared peak of galaxies at redshifts up to z ~ 4.5. The 
models use similar Monte Carlo fitting methods but dif- 
fer in a few ways, e.g. the BIO model is fit using data 
(predominantly number counts) from a very wide range 
of IR wavelengths (15 |J.m to 1.1mm) while MIO uses only 
data that constrains the evolution of the FIR peak and 
was not fit to any observations with A < 70 |xm. BIO 
also divides galaxies into two distinct populations based 
on luminosity and attempts to account for the strong 
lensing of high-redshift galaxies; MIO does neither of the 
above. 

2.4. Radio Galaxies 

Synchrotron and, to a lesser extent free-free emis- 
sion dominates the SEDs of radio galaxies at rest-frame 
A 1 mm. Thus radio galaxies become an increas- 
ingly important contribution to the CSB at wavelengths 
greater than ~ 1.5 mm {v ^ 200 GHz). Their number 

counts are relatively shallow (e.g.. lde Zotti et "aI]l2010D . 
meaning that their contribution to the power spectrum 
is dominated by the brighter sources, resulting in pri- 
marily Poisson noise, with the clustered term being sub- 
dominant. 

While radio sources are a significant source of power 
at 2030 and 1380 |J.m, they do not feature prominently 
in the cross-frequency correlation of ACT and BLAST 
maps. They only contribute to the uncertainties in the 
cross-power spectrum, so we do not include them in our 
models for the cross spectra. 
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Fig. 1. — Four of the maps used in this analysis: the BLAST maps at 350 and 500 ]xm (top) and the ACT data at 1380 and 2030 [om for 
the same region {bottom). Long wavelength modes in the ACT maps have been removed using the high-pass filter described in eqn. l|6)l. 
All maps are multiplied by a taper as discussed in the text. 



2.5. SZ 



Th e Sunyaev-Zeldovich effect (jSunvaev fc ZeldovichI 
|1972[ ) is the distortion of the microwave background due 
to interaction of CMB photons with free electrons in clus- 
ters, and consists of two main physical mechanisms, re- 
ferred to as "thermal " and "kinetic" . The thermal term 
(tSZ) is the inverse Compton scattering of CMB photons 
as they pass through the intra-clustcr medium in galaxy 
clusters. The result is that the intensity of the CMB 
spectrum longward of 1380 |J.m (218 GHz) decreases (i.e., 
a decrement), while that shortward of 1380 |J.m increases 
(i.e., an increment), and 1380 |J.m is the null for the non- 
relativistic case. The kinetic term (kSZ) is the Doppler 
shift of scattered CMB photons by the bulk motion of 
galaxy clusters. The strength of the signal is propor- 
tional to the product of the free electron density and line 
of sight velocity. 

2.6. Cirrus 

On large angular scales, a significant source of fluctu- 
ation power is emission from Galactic cirrus. Although 
the SEP field is among the least contaminated by cirrus 
in the sky (/moan = 1-16 MJy sr~^), contributions from 
Galactic cirrus must still be accounted for. 

The power spectrum of Galactic cirrus has been shown 
in many studies to exhibit power-law behavior. Its ampli- 
tude varies over th e sky, but its slope is always between 
2.6 and —3 fe.g.. iGautier et aI1ll992t iBoulanacr ct al.l 



Il996t iMiville-Dcschcncs et al.ll2007t iBracco ct al.ll2010D 

In the FIR/submm bands, the SED of Galactic cirrus 
is well described by a modified blackbody (Equation |3]), 
with T = 17.5 K and /3 ~ 1.9, peaking at A ~ 150 M,m 



(V09. [Bracco et ani2010f ). As a result, bands closest to 
the peak (in our case 250 |^m) are most susceptible to 
contamination. 

As o ne moves far from the SED pcak. llTnkbciner et al.l 
(|1999[ ) show that the modified blackbody approximation 
breaks down, and a multi-component fit is a much bet- 
ter description of the data. Therefore, for this analysis 
we measure the power spectr um at 100 M.m for the S EP 
region, and adopt model 8 of lFinkbeincr et ahl (|1999[ ) to 
estimate the amplitude of the power in our bands. 

3. INSTRUMENTS AND OBSERVATIONS 

Below we describe the ACT and BLAST data used for 
the cross-correlation analysis, as well as the data used 
for estimating the cirrus power spectrum. 

3.1. ACT 

ACT is a 6-m cter off-axis Gregorian telescope 
(| Fowler et al.ll2007[ ) situated at an elevation of 5190 me- 
ters on Cerro Toco in the Atacama desert in north- 
ern Chile. ACT has three frequency bands centered 
at 148GHz (2.0mm), 218GHz (1.4mm) and 277GHz 
(1.1 mm) with angular resolutions of roughly I'A, l.'O and 
0.'9, respectively. The high altitude site in the arid desert 
is excellent for mm observations due to low prccipitable 
water vapor and stability of the atmosphere. The trop- 
ical location of ACT permits observations on both the 
northern and southern celestial hemispher es. Further 
detail s on the instrument are presented in iSwetz et al.l 
(|2010H . lFowlcr et al.l (|2010[ ) and references thereii{3. The 



ACT Collaboration papers are archived 



at 



|http: //www .physics .princeton. edu/act/ I 
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ACT maps used in this paper are made from the 2008 ob- 
serving season data (148 GHz and 218 GHz, or 2030 and 
13 80 M.m, respec t ively) and are identical to th e maps used 
in lHaiian etall (|2010f ) and lDas et aP (|2010[ ). The beam 
full widths at half maxima (FWHM) are I A' and 1.0' a t 
148 GHz and 218 GHz, respectively (|Hincks et al.ll2010[ ). 
The maps have mean la sensitivities which vary slightly 
across the maps, ranging from 2.4 — 3.5 mJy beam"^ 
(median « 2.7 mJy beam^^), and 3.2 — 5.4 mJy beam~^ 
(median « 3 7 mJy beam'^) , at 148 GHz and 218 GHz, 
respectively ()Das et al.ll2010D . The map projection used 
is cylindrical equal area (CEA) with square pixels, 0.5' on 
a side. The ACT data-set is divided into four equal sub- 
sets in time, such that the four independent maps gen- 
erated from these subsets cover the same area and have 
si milar depths. We ca ll these "sub- maps" . As described 
m iHai lan cFaLl (|2010l) . the ACT maps are directly cali- 
brated to WMAP. This results in a 2% fractional temper- 
ature uncertainty for the 148 GHz maps. The calibration 
error for the 218 GHz maps is 7%. 

Because the ACT maps have poorly measured modes 
on the largest angular scales, we filter them using a high- 
pass filter Fc{£) in Fourier space. The high-pass filter is 
a smooth sine-squared function in Fourier space given by 

F,ii) = sin^ xi£)e{e - 4,in)e(Cax -e) + e{£ - w), 

(6) 

where x{£) = - £min)/(4iax - ^min) and 6 is 

the Heaviside function. We choose ^,„in = 100 and 
^max = 500. Moreover, the large-scale CMB in the 
ACT maps acts as noise in cross-correlations with the 
BLAST maps, since the CMB is absent in the latter. 
If not Itered, the large angular scale and CMB noise 
terms contaminate the real-space cross-frequency corre- 
lations described in Section [331 Therefore we use a filter 
with ^max = 2200 when dealing with real-space cross- 
frequency correlations. The analyzed power spectra are 
corrected for this filter as well as for the effects of the 
beam and pixel window functions. 

3.2. BLAST 

BLAST flew for 11 days from Antarctica in Decem- 
ber of 2006. Operating above most of the atmosphere, 
BLAST observed in bands which are difficult or not pos- 
sible to observe from t he ground. As a p athfinder for 
the SPIRE instrument (jGriffin et al.ll2003f) . it made ob- 
servations at 250, 350 and 500 \m\, of a number of tar- 
gets, both Galactic and extragalactic. Its 1.8-m under- 
illuminated primary resulted in beams with FWHM of 
36, 45 and 60 arcs econds. For a detai led d escription of 
the instrument see lPascale et al.l ()2008D and lTruch et al.l 

Among the fields BLAST observed is an 8.6 deg^ 
rectangle near the South Ecliptic Pole (SEP) — cho- 
sen because it is a relatively low-cirrus window through 
the Galaxy (see § 12. 6|) — whose corners lie at [ 
(5h02'",-52°50'); (4h57"^,-51°35'); (4^25"", -54°19'); 
(4h30"\-55°41')] (see Fig. [J). Further studies of the 
SEP field, includ uig BLAST ca t alogs and 2 4 M-m maps, 
can b e found in IValiante et al.l ()2010f ) and iScott et al.l 
(j2010( ). The maps have mean la sensitivities of 36.7, 

^ BLAST Results and P ublications can be found at 
|http : //blastexperlment . inf o"7] 



27.2, and 19.1 mJy beam^i at 250, 350, and 500 nm, re- 
spectively. Furthermore, confusion noise, due to multiple 
point sources occupying a single beam element, is esti- 
mated to be 7.6, 6.0 and 4.4 mJy beam~^ at 250, 350, 
and 500 |xm, respectively. The la uncertainty on the ab- 
solute calibration is accurate to 9.5, 8.7, and 9.2% at 250, 
350, and 500 |a.m, respectively. 

The BLAST time-ordered data (TODs) are divided 
into four sets — covering the same region of the sky 
to the same depth — from which we make four unique 
sub-maps. The number of subsets is chosen to maximize 
the number of maps that can be made while maintaining 
uniformity in hits and providing as much cross-linking 
as possible. The maps are made with th e iterative map- 
maker, SANEPIC (jPatanchon et al.l [20091) . resulting in a 
transfer function of unity on the scales of interest. These 
sub-maps are unique to this study and are publicly avail- 
able at http://blastexperiment.info/results.php. 

Due to poor cross-linking, however, large scale noise, 
resembling waves in the map, is present. This noise is 
easily dealt with by filtering in Fourier space, as de- 
scribed in § [4] Maps are made in tangent-plane projec- 
tion (TAN), with 10" pixels. In order to cross-correlate 
with ACT, BLAST maps are re-binned to ACT resolu- 
tion and reprojcctcd to cylindrical equal-area projection 
(CEA) using Montagel^ We confirm the alignment by 
analyzing the real-space cross-frequency correlation, de- 
scribed in detail in § 13.41 

3.3. IRIS 

To estimate the contribution from Galactic 
cirrus, we use three IRIS (reprocessed IRAS: 
iMiviUe-Deschenes fc Lagach^ [2001) HCO1N0 maps 
at 100 M.m. These maps are consistent with the 
iFinkb emer "cnH dHal EDS) ma p s used for estimating 
the Galactic cirrus in iDas et all (|2010f ). but they are 
at a higher resolution. Since we are most interested 
in large-scale modes, the power spectrum is measured 
for a 30 deg^ field surrounding the SEP (see Fig. [2]). 
The three HCONs are from the same region of the sky, 
with independent noise properties. The power spectrum 
of the Galactic cirrus is computed from the average 
cross-spectrum of these three maps. As a last step 
we correct the power spectra for a window function 
corresponding to IRAS's 4.3' beam. 

3.4. Comparing Data sets: Testing Alignment with 
Real-Space Cross-Correlations 

We test that the maps are properly aligned by inspect- 
ing their real-space cross-frequency correlations. This 
is done by inverse Fourier transforming the 2D cross- 
frequency correlation of the Fourier components of the 
maps: 

= Y.a{e)b*{£)F^{e)cxp{i£-^), (7) 

e 

where ^ is a vector in Fourier space, a{£) and b{£) are the 
ACT and BLAST maps in Fourier space, respectively. 



^ http: //montage . ipac . caltech.edu/ 

HCON refers to each individual observation at three different 
epochs. For more information, a nd publicly available maps, see 
|http : //www. clta.utoronto . ca/-ma md/IRIS/ 
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Fig. 2. — The full-sky IRIS map at 100 yxm scaled logarithmically. The red outline represents the BLAST field in CEA projection which 
is used for cross-frequency correlating with the ACT maps. The region clearly has low dust contamination. The bigger square area outlied 
by a green line shows the 30 deg^ region used for estimating the IRIS power spectrum. This map is filtered with the high-pass filter of 
Equation |6] 



and Fc{i) denotes high-pass filtering of Equation [S] in 
Fourier space. The resulting measured real-space cross- 
frequency correlation function M encodes the celestial 
correlation function between the bands, {C{x)) as: 



Maxfc(x) = C(x) ® B(x) + 7l(x), 



(8) 



where (g) denotes a convolution in real space, i?(x) is the 
effective beam between the two maps, and n(x) is the 
noise. Perfectly aligned maps would result in a 2D cross- 
frequency correlation function whose peak lies at x = 0. 
The shape of Maxb(x) depends on the correlation length 
of the field, t he hig h-pass filtering and the beams (See 
iHaiian "eFall (|20T0[ ) for a comparison between measure- 
ments of a real-space cross-freqency correlation function 
and simulations for a related example). As illustrated 
in Fig. [3l we measure a cross-correlation with a peak at 
x = (zero lag), indicating that the maps are correlated 
and properly aligned. 

4. MEASURING AUTO- AND CROSS-FREQUENCY 
POWER SPECTRA 

We examine three distinct types of power spectra: the 
BLAST auto-band spectra; the BLAST cross-band spec- 
tra; and the ACT/BLAST cross-frequency spectra. We 
also estimate the auto-band spectrum for the 100 |J.m 
IRIS maps. Details of our power spectrum method are 
given below. 

4.1. Power Spectrum Method 

We use three distinct types of power spectra: the 
BLAST auto-band spectra; the BLAST cross-band spec- 
tra; and the ACT/BLAST cross-band spectra. 

The maps can be represented as 



where ATsky(£c) is the sky temperature signal, N is 
the noise, B is the instrument beam and we use ® to 
represent a convolution in real-space. Both the ACT 
and BLAST maps are made with unbiased iterative 
map-makers, whose transfer functions are approximately 
unity on the angular scales of interest in this study, and 
can thus be safely neglected. 
All power spectra, both auto- and cross-frequency 




AT{x) = AT,ky(a;) (g) B{x) + N{x), 



(9) 
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Fig. 3. — The radial profile of the 2D cross-frequency correlation 
function between the ACT 1380 \mi band and the BLAST 500 yim 
is plotted with arbitrary normalization, with an image of the 2D 
function inset. The function has a clear peak at zero lage. This 
shows the two data sets are aligned and there is a correlated signal 
in the maps. Large scale fluctuations (at 9 > 500") are caused by 
the atmospheric noise and the CMB. The maps arc filtered by the 
high-pass filter of Equation [6] with ^niin = 100 and ^max = 2200 to 
remove longer wavelength noise. 
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Fig. 4. — The estimated 2D noise power spectrum for (a) 250 /^m and (b) 1380 fim. The noise model is computed from the difference of 
the auto- and cross-power spectra. Noise due to scan-synchronous signals and other large-scale correlations, which contaminates the signal 
in the central region, is down-weighted with this weight map in Fourier space, and then further filtered with a vertical mask spanning 
1^x1 < 500. The weights are normalized by the maximum value of the noise power spectrum and are plotted in a log-scale. 



spectra, are computed using cross-correlations of sub- 
maps (described in § 13.11 and § 13. 2p , the advantage being 
that the noise between sub-maps is uncorrelated and thus 
averages to zero in the cross-spectra. A cross-power spec- 
trum computed this way provides an unbiased estimator 
of the true underlying power spectrum. The power spec- 
trum methods used in this paper closely foll ow those used 
for cross-corr elating ACT and WMAP in Ha] ian et all 
(|20T?I . 

BLAST X BLAST power spectra are computed using 
the average of the six cross-correlations between the four 
BLAST sub-maps (in each band) , such that 

-, l</3<4 

^^=6 ^ (10) 

q,/3;q</3 

where a and /3 index the four sub-maps. The reason 
for six cross-correlations, rather than nine, is that cross- 
frequency cross-correlations of sub-maps made from the 
same scans are not used, in order to avoid introducing 
correlated noise or other systematic effects. 

Since the noise in all ACT and BLAST sub-maps is 
uncorrelated, we co-add the sub-maps for each frequency 
before computing the ACT x BLAST power spectra. The 
ACT X BLAST power spectra are computed from these 
maps. This is identical to cross-correlating all ACT sub- 
maps with all BLAST sub-maps and averaging them. 

Several components contribute to the cross-spectrum 
uncertainties. An analytic approach to computing the 
uncertainties is described in Appendix [BJ 

4.2. Weights and Masks 

We use a flat-sky approximation for computing the 
power spectra, which are measured from Fourier trans- 
forms of the maps. Since the maps are 2D, the power 
spectra that we obt ain are 2D. We ado pt several tech- 
niques developed in iHaiian et all (j2010D in order to ISO- 
late and remove or down-weight instrumental and sys- 
tematic noise. This is done in two stages, in real space 
and in Fourier space, before reducing the 2D spectrum 
into the familiar ID angular power spectrum binned in 
radius (£). The ACT x BLAST cross-spectra are limited 
by the area of the BLAST maps, which are 1.5 x 5.7° 



rectangles 8.6 deg^), rotated by approximately 30°, 
and with noisy edges (see Fig. [IJ . To apodize the sharp 
edges in the maps, we use the f irst Slepian tape r (prop- 
erties of which are described in lDas et al.ll2009l ) defined 
on the BLAST region. The gradual fall-off of this ta- 
per at the edges reduces the mode-mode coupling in the 
measured power spectrum. Any residual mode coupling 
caused by this weighting is corrected in the end by de- 
convolving the window function (i.e. the power spectrum 
of the taper function) from the me asured power spec - 
tru m using the algori thm described in lHivon" "crall(|200l 
and iDas et all ()2009D . Large-scale noisy modes are best 
treated in Fourier space. The statistical isotropy of the 
Universe leads to an isotropic 2D power spectrum from 
extragalactic sources, on average. Anisotropic power in 
2D Fourier space is caused by noise and is optimally dealt 
with using inverse noise weighting. This is done by di- 
viding the 2D spectra by our best estimate of the 2D 
noise power spectrum for each map. At each frequency 
band, the noise model is computed from the difference 
between the average 2D auto- and cross-spe c tra fo r each 
pair of maps as described in IHaiian et al.l (j201Clf ) . For 
every cross-correlation, noise weights are computed from 
the inverse of the square root of the product of the two 
noise power spectra corresponding to the two frequency 
bands. The weights are whitened by dividing by their 
angular averaged value with a fine binning. Using simu- 
lations we confirm that this weighting does not bias the 
signal power spectrum. Anisotropic, large-scale noise is 
evident at the center of the 2D power spectrum, as shown 
in Fig. |4l Noise-weighting down-weights the noisy ver- 
tical stripe that passes through the origin, which is pre- 
dominantly due t o large-scale unco nstrained modes in 
the map (see also lFowler et al.ll2010D . To further ensure 
that our results are not contaminated by this stripe, we 
mask a vertical band spanning \ix\ < 500 before aver- 
aging the 2D power spectra in annuli. Our results are 
invariant under further widening of this mask. In order 
to be consistent with the V09 analysis, we make a mask 
to remove all point sources that have a flux greater than 
0.5 Jy in the BLAST 250 iJ.m map (six sources). We use 
that mask for computing BLAST power spectra only. For 
cross-correlations of ACT and BLAST, we instead mask 
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Fig. 5. — Power spectrum of the 30 deg^ region surrounding the SEP from IRIS 100 ^im maps. A three component fit (left panel), 
including a clustering term (dashed yellow line) provides a better description of the measured power spectrum than fitting just Poisson and 
cirrus alone (right panel), but the difference is not statistically significant. 



just the brightest source in BLAST maps, which happens 
to be a local spiral galaxy. Masking more point sources 
does not have an effect on the cross-spectra. We a lso 
mask the known radio galaxies (jMarriage et al.ir2010[ ) in 
the ACT 2030 |J.m map to reduce the uncertainties in the 
cross-spectra (see Appendix [B|) . 

4.3. Estimating Galactic Cirrus Emission 

IRIS maps at 100 |J.m contain three potential sources 
of power - diffuse Galactic emission (or cirrus), as well 
as the Poisson and clustered terms of the DSFG power 
spectra. Cirrus dominates the power spectrum on angu- 
lar scales £ ^ 800, but varies depending on the observed 

patch of sky. The Poisson level is highly sensitive to 
the adopted flux cut. Thus, in order to detect the sig- 
nal from clustering, both the cirrus and Poisson noise 
must be sufficiently low. This is achieved by observing 
in a clean patch of sky, and cutting bright point sources. 
We realize these criteria by considering only the SEP 
(-?^mcan = 1-16 MJy sr~^), and by cutting sources with 

5' riit > 1 Jy- 

iMiville-Deschenes et al.l ()2007[ ) show that the power 
spectrum of Galactic cirrus can be approximated by a 
power-law, 

Pci„us(fce) = Po (^^) , (11) 

where kg is the angular wavenumber in inverse arcmin- 
utes, and Pq is the power spectrum value at fcg = 0.01 
arcmin" 

We are only concerned with the modes that affect our 
measurement, and since the cirrus power spectrum falls 
steeply with increasing kg, we focus our attention on the 
larger-scale modes. In order to probe these modes with 
maximum resolution in Fourier space, we measure the 



cirrus component of a ^ 30 deg^ region of the IRIS data 
surrounding the SEP field as indicated in Fig. [2] We fil- 
ter the IRIS maps using the high-pass filter of Equation|ni 
The power spectrum is computed from the mean of the 
three cross-spectra from the th ree HCON maps (using 
one t aper at resolution 1; see iDas. Haiian. fc Spergell 
[2001 for details). 

We attempt to fit the data in two ways, where in both 
cases we fix the slope of the cirrus power spe ctrum at a — 
—2.7 (adopting the properties of region 5 of lBracco et al.l 
[20T0I whose mean flux density most resembles the SEP). 
The first is a two-parameter fit, where the free variables 
are the Poisson level and the amplitude of the cirrus 
power. For this we find Pq = (0.47±0.06) x 10^ Jy^ sr"! 
and ~ 21.9 (dof = 29). The measured power spec- 
trum is shown in Fig. [Sj The power spectrum uncer- 
tainties are calculated in a manner analogous to that 
described in Appendix [BJ 

The second fitting procedure uses a three-parameter fit 
in which the free variables are the Poisson level, cirrus 
amplitude and clustering amplitude of the DSFGs. The 
shape of the clustering component is simply that of a 
linear dark matter spectrum. In this case we find cirrus 
values Pq = (0.19 ±0.15) x 10^ Jy^ sr~^ and a clustering 
amplitude of - 720 Jy^ sr~^ at £ = 3000, with 
18.1 (dof = 28). 

The two approaches estimate consistent Poisson lev- 
els, but the fit with a clustered component appears to 
describe the data better than without, with A^^ = 3.8. 
While not yet stati stically significant, fu ture studies with 
PACS at 100 i^m (jPoglitsch et al.ll2010D should be able 
to measure the clustered signal to high significance. 
The cirrus power spectrum is assumed to continue to 
smaller angular scales, and is estimated at the ACT 
and BLAST bands using the average dust emission color 
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Fig. 6.— BLAST x BLAST (250-500 \mi) and ACT x BLAST (1380-2030 \xm) power spectra in P{kg) with Ict errorbars. Squares and 
crosses are the data before and after cirrus removal. Red exes and red horizontal dotted lines are published power spectra and Poisson 
noise levels from V09. The dotted blue lines which are horizontal and which arc falling with kg, arc the best-fit Poisson and clustering 
terms, respectively. The departure from Poisson is the evidence for clustering of DSFGs. Note that the vertical scale is different for each 
panel. The error bars are described in Appendix [B] 



(-^(sub)mm/-^ioo)^; which is estimated using model 8 of 
iFinkbeiner eFall (|T999l ). 

5. POWER SPECTRUM RESULTS 

The BLAST auto-band and cross-band power spectra 
and BLAST x ACT cross-frequency power spectra are 
shown in Fig. [51 Raw data are shown as squares, while 
cirrus subtracted points are shown as crosses with error 
bars. The Galactic cirr us s pectra, interpolated to our 
bands as described in § 12.61 are shown as dashed lines 
in the bottom left corner of each panel (when strong 
enough to appear at all) . Cirrus appears to have a nearly 
negligible effect on the power in most bands, with only 
a marginal contribution in the 250 |J.m auto-spectrum. 
Note, the cirrus contribution in V09 to the BLAST 
bands was extrapolated from 100 [a.m incorrectly; how- 
ever, properly accounting for cirrus ultimately has little 
impact on the final result. The cirrus-corrected data are 
given in Table [5] We describe the models and the fits to 
these data in § [HI 

The figure shows a clear cross-correlation between 
ACT and BLAST. There is both a significantly corre- 
lated Poisson term (horizontal line) and a clear cluster- 
ing term (rising to low kg). This is the main result of this 
paper: that the unresolved BLAST background made up 
of DSFGs is intimately related to the ACT unresolved 
background. The signal is clearest in the ACT 1380 |J.m 
correlation with BLAST 500 and 350 \im., and less sig- 



nificant in the ACT 1380 and 2030 [a.m correlation with 
BLAST 250 iJ.m. Additionally, the figure confirms the 
V09 BLAST power spectrum analysis, and extends it to 
include the cross-frequency correlation between BLAST 
bands. 

Not shown in Fig. [5] arc predictions for the cross- 
correlation of the SZ increment and decrement, nor that 
of predictions for the cross-correlations of the SZ decre- 
ment and DSFGs. Both of these signals would appear 
as anti-correlatcd at the ACT 2030 |J.m band, and would 
act to decrease the total sky s i gnal. The former, us- 
ing templates of lBattaglia et al.l (|2010D . was predicted to 
be negligibly small; and while at some level the latter 
should exist, we have not yet identified a clear signature 
(which should appear only in the cross-correlations with 
the ACT 2030 ^m band). 

6. LINEAR CLUSTERING MODEL 

In this section we estimate the DSFG Poisson power 
levels, and fit the clustered component using a simple lin - 
ear model similar to that of V09 and iHall et al] ()2010[ ) . 
We assume that the clustered component of the DSFG 
power spectrum, PdsfGi is related to the linear dark mat- 
ter power spectrum, Pdm, through a single bias param- 
eter b{z): 

PDSFoik, z) = b{z fPuM{k, z), (12) 
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Fig. 7. — Redshift distributions of intensity, dl/dz, at 250 and 
2030 \aa for ten arbitrarily chosen realizations of the BIO source 
model. In this model the emission peaks around 2 = 1 for the whole 
range of wavelengths covered by ACT and BLAST, but there are 
significant contrilautions to the IR fiux from 2 2 at mm wave- 
lengths. 




2 - 



12 3 4 

Redshift 

Fig. 8. — Redshift-dependent best-fit linear bias h{z) for three 
realizations of the BIO model with Icr error bounds estimated from 
our MCMC. Also shown are the best-fit single-value linear biases 
found in V09 for the BLAST 250, 350, and 500 |jm auto-spectra, 
plotted at median redshifts z = 1.61, 1.88, and 2.42, respectively. 



SO that the angular power spectrum, P{kg), of DSFGs 
can be written as 

^^(z)^(z), (13) 

(jBond at all 119911 : iTegmark et"all |2002[ ). where x{z) is 
the comoving distance, dVc/dz = x^dx/dz the comoving 
volume element, and dl/dz is the contribution to the 
intensity from sources at redshift z. 

We adopt dl/dz from the source model of BIO. Fig. [7] 
shows plots of dl/dz from the model at 250 and 2030 \im 
for ten randomly chosen realizations provided by the 
BIO distribution^. We adopt the concordance model, 
a flat ACDM cosmology with VLm = 0.27 4, 17a = 0.726 
gp = 70.5 km s^^ Mpc"\ and cts = 0.81 (jHinshaw et al.l 
[200911 . 

The linear dark matter power spectrum is calculated 
as PDM(fc) = Poik)D^{z)T^{k), where Po(fce) is the pri- 
mordial power spectrum, T{kg) is the matter transfer 
functi on with fitting function given in lEisenstein fc Hul 
()1998D . and D{z) is the linear density growth function. 
For simplicity we treat the magnitude of the Poisson com- 
ponent of each of the 12 power-spectra as a free param- 
eter rather than simultaneously using the BIO model to 
predict the Poisson level as well as the clustering power. 

We note that our model does not account for non- 
linear, one halo clustered power. Though likely present, 
the data are not sufficiently constraining given that the 
Poisson power dominates over or is degenerate with the 
one-halo term in the angular scales to which we are sen- 
sitive. 

6.1. Estimating the bias 

In principle the bias, 6, is a function of scale 
and redshift, as well as environmental factors such 
as the host halo mass for the DSFGs. Here we 
adopt a simple redshift-dependent bias of the form 
(iBond. Carr. fc Hoganl[r991l: lHui fc Parfrevll2008D : 

&(z) = l + (&o-l)^, (14) 

where D[z) is the linear growth function and 6o is an 
initial bias at some formation redshift, zg. This parame- 
terization assumes that DSFGs are members of a single 
population, which formed at the same epoch (zq) and 
under the same conditions. 

Our parameter space consists of the 12 Poisson levels 
plus &o and zq- However, just as the Poisson contribu- 
tion is a sum over the galaxy distribution, weighted by 
the square of the fluxes, eqn. |4l so is the average bias , 
though weig hted linearly by the flux fe.g.. lBondlll993[) . 
Thus, in any physical model for the star forming objects 
and their bias the two terms would be correlated, but 
in a model-dependent way, so here we just adopt the in- 
dependent bias for simplicity. Decoupling the Poisson 
level and the clustered component means the interpreta- 
tion of our derived b is not straightforward. We find that 
moderate changes in zq (in the range 6 < zg < 10) have 
virtually no effect on the quality of the fit or best-fit 



® |http: //www. las .u-psud. fr/lrgal axles/model .php 



Correlations in the CSB from ACT x BLAST 11 



250 fj,m 350 fj,m 500 /um 1380 /um 2030 ixm 




1000 , 5000 1000 5000 1000 5000 



Fig. 9.— BLAST x BLAST (250-500 ^m) and ACT x BLAST (1380-2030 \mi) power spectra in i^Cijl-K. Data, which have had Galactic 
cirrus power removed, are shown as blue crosses. Red squares are the same data after removal of the Poisson term, and after logarithmic 
binning, with log(A£) = 0.2, and represent the contribution to the total power spectrum from clustering. Pink exes are the clustered term 
data from V09. The blue dotted lines rising to larger £ are the best-fit Poisson terms, and the approximately horizontal blue dotted lines 
are the best-fit clustering terms, which are determined by the 2-dependent bias, as described in §(6] Also plotted are the phenomenological 
models of BIO and MIO, in green and brown (with shaded error regions), respectively. Poisson levels are calculated after truncating the 
counts at 300, 250, 170, 20, and 20 mjy at 250, 350, 500, 1380 and 2030 |jm, respectively. Error regions are calculated with Monte Carlos. 
Both the models agree at some effective wavelengths, but disagree at others, so that neither describes the data fully. The MIO model also 
somewhat over-predicts the CIB at BLAST wavelengths, which is consistent with the behavior of the model Poisson term here. Note that 
the vertical scale is different for each panel. The cirrus-corrected data here are given in Table [2] 



and also find that and zq are almost degener- 
ate, and so we fix zq = 8. We explore the remaining 
13-parameter space using a Markov Chain Monte Carlo 
(MCMC) method with uniform priors on each parame- 
ter. We fit the ACT X BLAST data in the range i = 950 
(ke = 0.04 arcmin-i). We find a best-fit 6o = 18.2^^;^ 
with = 107 for 101 degrees of freedom. This corre- 
sponds to h{z = 1) = 5.0l[Jj and 6(z = 2) = 6.8t;^j. 
This Poisson plus clustering model is preferred to the 
null case with no ACT x BLAST correlation at over 2ha. 
The best-fit clustering and Poisson levels are reported in 
Table [H and are plotted as blue dotted lines in Figs. IH] 
andO 

We also try fitting just the ACT x BLAST data with 
no clustering power. In this case we obtain a best-fit 
= 64.3 with 48 degrees of freedom. After adding 
linear clustering to the model, we find = 43.6 with 47 
degrees of freedom, corresponding to Ax^ of 20.7 (with 
one fewer degree of freedom) , so that the model including 
clustering is preferred to one with only Poisson and cirrus 
at greater than 4(T. Additionally, the Poisson levels are 
lower when clustering is included. 

Lastly, we try fitting a single-value bias, independent 
of redshift, to the entire range of power spectra. We find 



a best-fit single-value bias of 5.0±0.4 with x^ = 110.8 for 
101 degrees of freedom; worse than a redshift-dependent 
bias by Ax^ = 3.8. The single-value bias thus provides a 
good fit to the ACT x BLAST data, however when we in- 
clude the measured 2 030 jam (ARl) and 13 80 ^m (AR2) 
clustered power from lDunklev et"an ()2010D in the likeli- 
hood calculation, wc find the redshift-dependent bias is 
preferred at ~ 2(7. The ACT clustering measurements 
are reproduced well with h{z) (see Fig. [TOl right panel) 
but under-predicted with the single-value bias. 

Fig.[S]shows hiz) calculated using Equation[T3]for three 
realizations of the BIO model which span the entire range 
of results. The bias appears high compare d to the linear 
bias estimates in V09 (who adopted the iLagache et al.l 
[200l model), although given the spread in h[z) from 
different realizations of the BIO model, the measurements 
are not inconsistent. 

Our choice of source model and bias parameteriza- 
tion is likely affecting h{£). The BIO model contains 
two distinct classes of IR sources, "normal" and "star- 
burst," with substantially different luminosities, which 
is not accounted for in our 6(z) parameterization. Also, 
the BIO model does not match observational constraints 
equally well across the whole wavelength range probed 
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by ACT and BLAST; for instance it under-predicts Pois- 
son power and number counts compared to BLAST and 
SPT measurements a t 500 and 1360 |i.m (see § 5.6 of 
iBethermin et aLll2010L as well as Fig. |9] of this paper). 
An under-prcdiction of dl/dz would result in a higher 
bias to compensate. 

7. DISCUSSION 

We compare the models of BIO and MIO (described in 
§ 12.31 ) to our results. Poisson levels are calculated from 
the number counts using Equation |4l where the values 
of S'cut are chosen to mimic those used in the analysis, 
i.e., Scut = 500, 250, 170 mJy at 250, 350, and 500 nm, 
respectively. 

Model predictions of BIO and MIO are shown in Fig. |9] 
as green and brown lines with shaded error regions, re- 
spectively. Both models agree at some effective wave- 
lengths, but disagree at others, so that neither appears 
to fully describe the data. The MIO model also somewhat 
over-predicts the CIB at BLAST wavelengths, which is 
consistent with the behavior we see for the model Poisson 
term. On the other hand, as already mentioned, the BIO 
model under-predicts Poisson power and number counts 
compared to BLAST and SPT measurements at 500 and 
1360 |J.m, which again is consistent with the behavior of 
the Poisson term of the model in Fig. [9l 

The Poisson and clustered power amplitudes are plot- 
ted as a funct ion of effective wavelength, defined as 
Acff = \/AiA2, in Fig. [101 Also included are mea- 
surements made by the following experiments: AKARI 
at 90 i^m (iMatsuura et al.l I2010D : Spitzer at 160 |i.m 
Lagache ct al.ll2007D : BLAST at 250, 350, and 500 M-m 



Viero et al.l 120091): ACT at 1380, 1673, and 2030 nm 



Dunklev et al.ll2010D : SPT at 1363, 1629, and 1947 ^m 
Hall ct all 120101 ): and the FIRAS modified blackbody 



{T ~ 18.5, /3 = 0.64), which is shown as a dotted line. 

The degree of correlation between widely spaced wave- 
lengths is of interest both in determining the redshift 
distribution of sources, and for modeling the IR source 
power as a CMB contaminant. To assess the correla- 
tion, the geometric means at each effective cross-band 



wavelength, defined as y • G^' , arc shown as a down- 
ward pointing arrows. Since we do not measure at 
A = 1380 and 2030 \xm, we rely on measurements by 
iDunklev et al.l ()2010D for those bands when calculating 
the geometric means. The ratios of the measurements to 
the geometric means, C^^ /{C^ ■ )^/^, then represent 
the levels of cross-correlation between bands. These are 
shown in Fig. [TT] for the Poisson power as a function both 
of effective wavelength and of distance between bands. 
Correlation is seen between all the frequencies, and does 
not fall significantly as a function of increased band sep- 
aration, suggesting a tight redshift distribution for the 
overlapping popula tion. This behavi or is consistent with 
the fin dings of e.g.. iHall et al.l ()2010[ ) and IDunklev et all 
()2010D : that the 1000 - 2000 |^m bands are correlated at 
close to the 100% level, and extends the range of wave- 
lengths probed. 

8. CONCLUSION 

We present measurements of the auto- and cross- 
frequency correlations of BLAST (250, 350 and 500 |^m) 



and ACT (1380 and 2030 ^im) maps. We find signifi- 
cant levels of correlation between the two sets of maps, 
indicating that the same DSFGs that make up the un- 
resolved fluctuations in BLAST maps are also present in 
ACT ma ps. Furthermore, we confirm previous BLAST 
analyses (jViero et al.ll2009l ) for a different field and with 
an independent pipeline, and extend the analysis by in- 
cluding BLAST X BLAST cross-frequency correlations. 

We fit Poisson and clustered terms at each effective 
wavelength simultaneously, which we achieve by ad opt- 
ing a model for the sources ([Bethermin et al.]l2010( ). as- 
suming a parameterized form for the z-dependent bias 
and using an MCMC to minimize the x^- Using this 
model we detect a clustered signal at Aa, in addition to 
a Poisson component. The best- fit bias is one that in- 
creases sharpl y with redsh i ft, an d is consistent with what 
was found bv I Viero etlTI (|2009D . 

We compar e phc nomeno logical mode ls by 
IBethermin et all ()2010D and iMarsden et"an (|2010f ) 
to the data and find rough agreement at numerous 
effective wavelengths. But we also find that neither 
model quite reproduces the data faithfully. Thus, 
we expect this measurement and others like it will 
ultimately provide powerful constraints for the redshift 
distribution and SEDs of future versions of the models. 

Though we find convincing evidence for correlated 
Poisson and clustered power from DSFGs, the levels of 
precision needed to robustly remove these signals from 
CMB power spectra demand better measurements still. 
This is particularly true of the clustering term, whose 
contribution to the power spectrum in PGi peaks at 
i ~ 800 - 1000, which is also the region in £-space typ- 
ically targeted in searches for the SZ power spectrum. 
Since the clustered term should scale independently of 
the Poisson term, the measurement becomes increasingly 
important to determine precisely. Future studies combin- 
ing Herschel/SVmE with ACT, SPT, and Planck wiU 
go a long way towards solidifying this much needed mea- 
surement. 
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Fig. 10. — Ci versus wavelength for observations and models. From left to right, the actual or effective wavelengths, X^ff = y/\iX2, 
(invrni) arc: 2030, 1673, 1380, 1007, 843, 831, 712, 695, 587, 500, 418, 354, 350, 296, 250, 160, 100, and 90. Best-fit Poisson (left panel) and 
clustered (at £ = 3000, right panel) Ci from measurements are shown as squares (auto-spectra) and diamonds (cross-spectra), respectively, 
and our measurement of IRIS galaxies arc shown as red crosses. Open circles represent the prediction for the clustered power at the ACT 
wavelengths from the best-fit, redshift-dcpendent bias model. Uncertainties are omitted for visual clarity, but are generally smaller than the 

size of the symbols due to the large dynamic range in . The geometric mean of the cross- band spectra, defined as ^ C^^ ■ C^^ , are shown 

as do wnward-pointing arro ws. Measurements fro m other experiments ar c: ACT (IDunklcv ct al.l2010l . yellow asterisks ); BLAST ll Viero et al.l 
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black triangle). The FIRAS modified blackbody (T = 18.5, /3 = 0.64) is plotted as a dotted line. As was seen in lHall et al.l 112010 . Fig. 5), 
FIRAS describes the data short of 500 |jm, but over-predicts the measurements at millimeter wavelengths. The ratio of the measurement 
(diamonds) to the geometric mean (downward-pointing arrows) represents the level of cross-correlation between bands. 
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TABLE 1 

Best-fit c'f°"'''°" and f7^i"stermgj.^ _ 3000-)^ including predictions for the clustered power at the three effective ACT bands. 
Predictions for the Poisson power at ACT bands are not provided as the Poisson terms are treated as free parameters 

when obtaininc; the best fit (see §[6ll. 
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TABLE 2 

Measured {/jK'^) from Fig. [9] (blue crosses). 
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Fig. 11. — Left panel; Cross-frequency correlation versus effective wavelength for Poisson. Cross-frequency correlation is defined as 
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APPENDIX 
A. UNIT CONVERSION 

The flux density unit of convention for infrared, (sub)millinieter, and radio astronomers is the Jansky, defined as: 

Jy=10-2«^— , (Al) 
m^ Hz 

and is obtained by integrating over the solid angle of the source. For extended sources, the surface brightness is 
described in Jy per unit solid angle, for example, Jy sr~^, (as adopted by BLAST), or Jy beam~^ (e-g-, SPIRE). 
Additionally, the power spectrum unit in this convention is given in Jy^ sr~^ 

CMB unit convention is to report a signal as STqmb', the deviation from the primordial 2.73 K blackbody. To convert 
from Jy sr^^ to STqmb in /^K, as a function of frequency: 

ST.= ('-^), (A2) 



ST 



SB^ 2kfkTcMBV a^^e^ 98.91Jy sr"! x^e^ 

(|Fixsenl [20091 ). Because the BLAST bandpasses have widths of ^ 30% (jPascale et al.l 120081 ) . and because the CMB 
blackbody at these wavelengths is particularly steep (falling exponentially on the Wien side of the 2.73 K blackbody), 
the integral of SB^/6T over the bands is weighted towards lower frequencies; an effect that becomes dramatically more 
pronounced at shorter wavelengths. Thus, the effective BLAST band centers in 5T arc ^ 264, 369, and 510 jam, leading 
to factors of conversion from nominal of 2.46, 1.75, and 1.13, respectively. 

Lastly, the CMB power spectrum is conventionally reported versus multipolc while in the (sub) millimeter the 
convention is to report it versus angular wavenumber, kg = 1/0, which is also known as a in the literature, and is 
typically expressed in arcmin"^. In the small- angle approximation the two are related hy I — 2'Kkg. 

B. POWER SPECTRUM UNCERTAINTIES 

The contents of each map used for cross-frequency correlations can be considered as a sum of two parts: one with a 
finite cross-correlation; and the other with vanishing cross-power spectrum. The former contributes to the signal in the 
cross-power spectrum, while the latter contributes to the uncertainties. Therefore three terms contribute to the power 
spectrum uncertainties: sample variance in the signal due to limited sky coverage; the noise; and a non-Gaussian term 
due to the Poisson distributed compact sources and galaxy clusters. The diagonal component of the ACT x BLAST 
cross-spectrum variance can be written as the sum of these terms, in order: 



nb\I Tib /sky 

where Nb^ estimates the average power spectrum of the noise; the superscripts and label the maps (1 for ACT, 
2 for BLAST); rtf, counts the number of Fourier modes measured in bin b (that is, the number of pixels falling in the 

appropriate annulus of Fourier space); /sky is the patch area divided by the full-sky solid angle, 47rsr; and Cb is the 
mean cross-spectrum. The last term, cTp, arises from the Poisson-distributed components in the maps (i.e., unresolved 
co mpact source s and c lust ers of galaxies ) and i s given by the non-Gaussian part of the four-point function as described 
in lFowler et al.l ()2010D and lHaiian ct al.l ()2010D . For purposes of the covariance calculation, we assume that the spatial 
distribution of these objects is uncorrelated. This term is constant with £. 
The noise terms in the ACT and BLAST maps arc given by 

7Vf=<, (B2) 

where C^'^ is the power spectrum of the radio galaxies in the ACT maps and and are the noise spectra in 
ACT and BLAST respectively. The noise ter ms, Nb, are domm ated by the atmospheric noise on large angular scales 
and by detector noise on the smallest scales (jPas et al.ll2010( ). The CMB is a major source of noise for this study 
out to ^ ^ 2500, especially for the 148 GHz data. Radio galaxies only contribute to the uncertainties through the 
fourth moment of the field. They do not bias the signal. The effect of the radio galaxies is stronger at 148 GHz and is 
negligible at 218 GHz. Therefore we mask the brightest radio galaxies in the 148 GHz map to reduce the uncertainty 
on the cross-power spectra. 

The unc ertainties on BLAS T x BLAST power spectra arc computed using a similar analytic estimate, given in 
Eqn.(9) of iFowler et"all (|2010f ) with riu, = 6 cross-spectra per map. 
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As a sanity check, we compare our analytic estimate of the error bars with the standard deviation of the power 
spectra computed from patches of the sky. We divide the data into four patches of equal area and with them compute 
four independent cross-power spectra. We use the variance of the measurements at each (. bin as a measure of the error 
on the power spectrum. This method agrees well with the analytic estimate of the errors; however, due to the small 
area of the sky used in this analysis, both analytic and patch- variance estimates of the error bars have uncertainties 
which arc limited by /sky Thus, we conservatively use the greater of the analytic and patch- variances as an estimate 
of the uncertainty of the power spectrum. We test the effect that this choice of error bars has on our results in section 

When fitting p arameters, we take the joint likelihood function to be diagonal as the off-diagonal elements are small 
(iDas et al.ll2010l ). 



